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We theoretically investigate quantum-mechanical dynamics of quasi-one-dimensional boson- 
fermion mixtures of atomic gases trapped in a toroidal potential, where effective inter-atomic inter- 
actions are tunable and affect the dynamics. We especially focus on effects of quantum statistics and 
many-body correlations beyond the Hartree-Fock (HF) mean-field approximation on the dynamics. 
In order to predict the dynamics, we utilize the numerical exact diagonalization method and also 
reproduce the calculation in the HF approximation for comparison. The toroidal gases originally 
have a rotational symmetry in the toroidal direction. We firstly prepare a deformed ground state 
as an initial state by adding a weak potential deformed in the toroidal direction, and then remove 
the potential to start the dynamics. In the dynamics, number densities of the deformed gases ex- 
hibit oscillations as demonstrated in the present paper. As a result, we find out that the bosons 
and fermions show quite different behaviors owing to quantum statistics. In particular, the bosons 
exhibit a low-frequency oscillation in the strong boson-boson attraction regime owing to the many- 
body correlations, and it can not be reproduced in the HF approximation. The oscillation of the 
fermions is strongly influenced by that of the bosons through the boson-fermion interaction as a 
forced oscillator. In addition, we also discuss a relationship between the low-frequency oscillation 
and restoration of the broken symmetry. 

PACS numbers: 67.85.-d, 67.85.Pq 



I. INTRODUCTION 

According to quantum statistics, particles are clas- 
sified into two types, i.e., bosons and fermions obey- 
ing the Bose-Einstein and Fermi-Dirac statistics, respec- 
tively. Bosons may occupy the same single-particle state 
and can produce the Bose-Einstein condensation, where 
macroscopic numbers of bosons occupy a single-particle 
state. On the other hand, fermions may not occupy the 
same single-particle state and can produce the Fermi de- 
generacy, where each fermion occupies a single-particle 
state up to the Fermi level. Quantum statistics plays an 
important role in quantum many-body physics not only 
for infinite matter systems but also for finite quantum 
systems. 

Recent developments of cooling and trapping tech- 
niques of atoms enable us to investigate quantum- 
statistical properties of ultra-cold atomic gases. In fact, 
atomic Bose-Einstein condensates [jj, Fermi degenerate 
gases Q , and boson-fermion mixtures are realized 

and studied in actual experiments. 

As a feature of the ultra-cold atomic gases, precise op- 
tic and electromagnetic techniques have produced various 
quantum systems designed at will. In particular, the Fes- 
hbach resonance technique gives tunable effective inter- 
atomic interactions from attractions to repulsions |13| . In 
addition, designed trap potentials have realized artificial 
crystal structures in the optical lattices and quasi-low- 
dimensional gases in the squeezed and toroidal poten- 
tials [3. 

Since bosons and fermions exhibit the quite differ- 
ent behaviors in the ultra-cold atomic gases, the boson- 
fermion mixtures present rich physics owing to the dif- 



ferent quantum statistics and boson-fermion interactions. 
In fact, as static properties, rich phase structures are pre- 
dicted in the atomic boson-fermion mixtures fl5| . In ad- 
dition, atomic fermions can produce the boson-fermion 
mixtures of atoms and molecules with quantum many- 
body correlations among atoms and molecules [HI, [l7[ . 
Furthermore, as dynamical properties, collective excita- 
tions [H, [lj| and collapse phenomena @, [2(J are studied 
theoretically and experimentally. 

In recent years, the ultra-cold atomic gases in the 
toroidal trap potential, called toroidal or ring-shaped 
gases, have been studied theoretically and experimen- 
tally. In fact, the toroidal bosons are experimentally 
demonstrated with the magnetic field 21]. the optical 
field l2p , the magneto-optical field [23[ , the microchip 
ring [24]], and the ring- lattice 1251 . In addition, theoretical 
works on the toroidal bosons |26l |27| and boson-fermion 
mixtures [28| are also published. 

As a feature of the toroidal gases, the gases can ex- 
hibit quasi-one-dimensional behaviors under a periodic 
boundary condition owing to geometry of the potential. 
Thus studies on the toroidal gases are closely connected 
to those on the one-dimensional systems. On another 
front, the toroidal gases are obviously a type of finite 
quantum many-body systems. In this sense, studies on 
the toroidal gases can also offer an important infrastruc- 
ture for the finite quantum many-body problem. Our 
aim of the present study is near the latter side. 

In the quantum many-body problem, dynamical prop- 
erties are fundamentally important. That is because the 
dynamics can give us significant information of quantum 
systems not only in theories but also in experiments. In 
fact, an indication of a spontaneous symmetry breaking 
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(or phase transition) in the thermodynamical limit can 
be observed in the corresponding dynamics in finite quan- 
tum systems as demonstrated in the present paper and 
also, e.g., in nuclear physics [29j . 

As an elementary knowledge of the quantum many- 
body physics, the mean-field approximation, which is 
conventionally used in studied on the cold atoms, gets 
worse near the phase transition in infinite matter sys- 
tems and also near the corresponding situation in finite 
quantum systems. Furthermore, in parallel, the mean- 
field approximation also grows worse as the interactions 
get stronger. In fact, it is theoretically confirmed for the 
toroidal bosons [26| . 

In the present paper, we theoretically study the dy- 
namics of boson-fermion mixtures of the toroidal gases 
as a finite quantum many-body problem linked to the 
actual experiments. In the dynamics, we especially focus 
on effects of quantum statistics and many-body correla- 
tions beyond the mean-field approximation. 

In order to predict the dynamics, we utilize the nu- 
merical exact diagonalization method and also reproduce 
the calculation in the Hartree-Fock (HF) mean-field ap- 
proximation for comparison. This procedure must be a 
conventional scheme to treat finite quantum systems and 
give exact solutions except for numeric errors and limi- 
tation of the configuration space on our computer mem- 
ories. 

As explained above, the present study is closely con- 
nected to the theories of the one-dimensional systems. 
The theories are developed in mathematical physics [30| . 
Especially those of the solvable models are studied well 
in the Bcthc ansatz method 
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o ne-dimensional 
bosons boson-fermion mixtures [33j, and so on. In 
addition, those of the soliton solutions are also studied 
in the mean- field approximation |34| . 

This paper is organized as follows. In Sec. [TTl we pro- 
pose a theoretical formulation used in the present paper. 
In Sec. lIIII we show static properties of the toroidal gases 
as a background knowledge to study the dynamics. In 
Sec. IIV1 we demonstrate the dynamics in some typical 
situations. In Sec. [V] we discuss the results obtained in 
the previous section. Section [VII is devoted to summary 
and perspective. 



II. FORMULATION 

Here we formulate a theoretical model for the toroidal 
gases in Sec. Ill Al and its quantum-mechanical dynamics 
in Sec. IIIB1 and also explain the calculational schemes 
used in the present paper in Sec. Ill CI In particular, the 
deformation parameters calculated in the next section are 
introduced in Sec. Ill Bl 

In order to simplify the following descriptions, we de- 
termine the system of unit as h = 1 for the reduced 
Planck constant without loss of generality. 



A. Toroidal gases 

The toroidal gases are the dilute atomic gases trapped 
in the toroidal potentials in vacuum, where the typical 
size of the gases is much greater than that of the atomic 
scale, or the Bohr radius. 

In the cooling techniques of atoms, the gases are cooled 
down near zero temperature and exhibit purely quantum- 
mechanical behaviors. 

First we explain a theoretical treatment of the dilute 
atomic gases, and then derive the quasi-one-dimensional 
effective Hamiltonian in a conventional manner. 



1. Theory of the dilute atomic gases 

The atoms of the dilute gases mostly behave as free 
particles except for extremely-short effective ranges of the 
inter-atomic interactions. In fact, the effective ranges are 
comparable with the atomic scale, and mean free paths 
of the atoms are much greater than the atomic scale in 
the dilute gases. 

In quantum mechanics, the many-body wave-function 
of the atoms exhibits singular behaviors near two-atom 
contact points owing to the extremely-short-range inter- 
actions, and the singularities themselves are mostly inde- 
pendent of gaseous behaviors because of the great differ- 
ence of the length scales. The singularities mostly agree 
with those in the inter-atomic scattering problem. 

In order to simplify descriptions for the gaseous behav- 
iors, we consider an asymptotic wave-function, instead of 
the real many-body wave-function, in the following man- 
ner. First we theoretically assume the asymptotic wave- 
function derived from the real wave-function by getting 
rid of the singularities. Second we renormalize the inter- 
atomic interactions with the corresponding low-energy 
scattering data, where information of the singularities 
are effectively included in the scattering data. 

The asymptotic wave-function approximately agrees 
with the wave-function of the non-interacting particles 
except for non-singular effects of the renormalizcd inter- 
atomic interactions. As a result, we can utilize the 
perturbation theory and mean-field approximation for 
the asymptotic wave-function unlike the case of the real 
wave- function. 

In the notations of the second quantization, we use the 
asymptotic field operators corresponding to the asymp- 
totic wave-function, instead of the real field operators, in 
a conventional manner 13511 . 



2. The quasi- one- dimensional effective Hamiltonian 

Because the asymptotic fields roughly agree with those 
of the non-interacting particles as explained above, we 
first consider the toroidal gases of the non-interacting 
particles, and then introduce the renormalized effective 
interactions in a conventional manner. 
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Let us consider the mixed gases of the single- 
component bosons and fermions in the toroidal potentials 
denoted by 



v a (r) = vi c \r c )+v^>(r z ) 



(*)/ 



(1) 



for a — b and f, where the subscripts b and f indicate the 
potentials for the bosons and fermions, respectively. The 
potentials v a (r) in Eq. (JXJ) are separated into the trans- 
verse and longitudinal potentials, Va\r c ) and Va (r z ), 
in the cylindrical coordinates, r = (r c ,9,r z ), and inde- 
pendent of the azimuthal angle 9. 

The transverse potential Va\r c ) in Eq. (P) has a deep 
valley centering at r c = R with the toroidal radius R, 
and can be described as 



Rf 



(2) 



with the transverse trap frequencies oJf and ojf near 
the center of the valley. Here we introduce the boson and 
fermion masses and m,f, respectively. Since the gases 
are mostly trapped in the valley, the detailed profiles of 
the transverse potentials hardly affect the following for- 
mulation. 

The longitudinal potential Va\r z ) in Eq. (JTJ) can be 
denoted by 

vi z \r z ) = \m a ^fr\ (3) 

(z) (z) 

with the longitudinal trap frequencies and ojf . 

The Hamiltonian for the non-interacting bosons and 
fermions is denoted by 



H a = H h + H { 
with the boson and fermion parts, 



Hn 



dr ipl(r)h a (r)il; a (r) 



(4) 



(5) 



for a = b and f, where we introduce the single-particle 
Hamiltonian 



h a (r) 



1/19 d Id 2 d 2 



2m a \r c dr c dr c r 2 d6 2 dr z 2 



-v a [r 



(6) 

The asymptotic field operators tpbir) and ip{( r ) f° r the 
bosons and fermions, respectively, obey the commutation 
and anti-commutation relations, 

[^(r),Vb(*)]_=0, [Mr)^l(s)}_ = 6(r-s), (7) 



[0f(r),^(*)] + = O, [Mr)^l(s)] + = 6(r-s), (8) 
and 

[lMr),V*(*)]_=0, [Mr),$(s)]_=0- (9) 



The eigen-states of the single-particle Hamiltonian 
h a {r) in Eq. © can be described as 

[K{r) - e mmK \ ^ a c L(rM z l(r z )^\0) = (10) 

with the eigen-energy e anmK , radial part 4>al K (r c ), Ion- 
gitudinal part <fiani(r z ) and azimuthal part 

(11) 



V2tt 

for the quantum numbers n (> 0), m (> 0), and n. These 
parts are orthonormal as 



dr c r c 



dr. 



&L(rc)]*<t> ( :L(rc) = 6 n , n ,, (12) 
^l(r z )Y^l(r z ) = 6 m , mr , (13) 



and 



d0 



4>W{e) 4> K :.'{e) = 6 K 



(14) 



By using complete sets of the eigen-states in Eq. (|10l) , 
we expand the field operators as 

Mr) = ^\r) + ^\r) (15) 

with the toroidal part 



E 

K — — 00 



&(rc)^(rz)^ e) (0)a a oo K (16) 



and residual part 



oo oo 



tf'W^EE E ^L(rc)^l(r z )^H0)a anmK - 

n—1 rn—1 k,— — oo 

(17) 

According to Eqs. the particle operators, ab nmK 

and af nmK , obey the commutation and anti-commutation 
relations, [a hnmK , al n , m , K ,]_ = <5„,„'<5m,m'£ K ,K' and so on. 

In the toroidal gases, although the interactions break 
the variable separation form in Eq. (ITU1) . the particles 
mostly occupy the n = m = states in the asymptotic 
field picture. That is because the excitation energies, 
£aWK - £ a ooK and e a oi K ~ e a0 o K , are much greater than 
both of the Fermi energy and absolute values of the in- 
teraction energies. 

At that time, we can focus on the limited configuration 
space described by ip^ (r) and (r), instead of the 
full configuration space, as an effective theory. This the- 
ory can be reduced to a type of the effective interaction 
theory pKjj ]. 

The toroidal parts ip\^\r) and ip\ P \r) in Eq. (ITBl can 



be rewritten as 



OO 

E ■ 

K — — OG 



i e \6)b K + 0(R- 3 ) 
(18) 
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and 



and 



K, — — 00 



with b K = aboo/c and c K = amo K . Here 



(<=)' 



1/4 



m a uj a 
ttR 2 



-m a u^ ) {r c -R) 2 /2 



and 



m a uj K a 



1/4 



-m a u4*VJ/2 



(19) 



(20) 



(21) 



In the limited configuration space, the ffb and fff in 
Eq. ([5]) are effectively reduced to 



and 



fff 



E 



K — — OO 



2m h R 2 



blb K + 0(R- 3 ) + £ b (22) 



E 



K, — — OO 



2m f i? 2 



4c K + 0(R~ 3 ) + St (23) 



with the constant energies £b and Sf. The constant ener- 
gies do not affect physics in the nonrelativistic systems. 

Because the higher terms of the R inverse, 0(R~ 3 ), 
must have negligible contributions to physics of the 
toroidal gases, we neglect them in the following formula- 
tion. At the same time, the quantum number k should 
be limited in |k| < n c with a large cut-off number k c in 
order to justify the approximations in Eqs. (TT51) and (111)1) . 

As a result, the effective Hamiltonians in Eqs. (|2"2")l and 
(l23l) indicate those of the quasi-one-dimensional gases of 
the non-interacting bosons and fcrmions. 

Next we consider the renormalized effective interac- 
tions -ffbb and ffbf in the limited configuration space, 
where ffbb and ffbf indicate the boson-boson and boson- 
fermion interactions, respectively. Note that the interac- 
tion between the fermions is neglected owing to the Pauli 
blocking effect. 

In a conventional manner, we utilize the renormalized 
contact-type pseudo-potentials for the low-energy s-wave 
scatterings, i.e., 

H hh = 9 -f fdr4 P) \r)4 P) \r)4 P \r)^ h P \r) (24) 



and 



ffi 



bf — ffbf 



dr ipl P) \r)4 P) \r)4 P) (r)ipl P) (r). (25) 



27r(mb + mi) 

gbi — abf 

m b mf 



(27) 



with the boson-boson and boson-fermion s-wave scatter- 
ing lengths abb and abf, respectively. The proportional- 
ity relations in Eqs. (f2"tj|) and (|27l) must be valid for the 
weakly-interacting gases as confirmed in the perturba- 
tion theory and mean-field approximation, or the Bethe- 
Goldstone theory, and should be verified for the strongly- 
interacting gases experimentally. 

It should be noted that the contact-type interactions in 
Eqs. ([24]) and (|25|) can not produce divergence in the two- 
body correlations because of the limited configuration 
space unlike the case of the full configuration space. The 
divergence in the full configuration space originally comes 
from the singularities owing to the extremely-short-range 
interactions and must be removed in the asymptotic field 
picture. 

The interaction parts ffbb and -ffbf in Eqs. (|24|) and 
(f2"5")l are reduced to 



ff 



bb 



and 



ff 



bf 



ffbb 
AttR 



ffbf 
2vri? 



E E E E 5 («i+«2),(k 3+K4 ) & k 1 & L^3^ 



ftl ft 2 K3 K4 



EEEE*( 

Ki K-2 ft 3 ft 4 



(28) 

2 ),(k3 + K4)^Ki C K2 Ck 3^ k <1 

(29) 



with the redefined coupling constants 



ffbb 



, (c) (z)~ 

m bV w b w b 9bb 



and 



ffbf ; 



(c) (z) (c) (z)~ 

m b TO f A/ w b 'w h 'w\ 'u\ 'gbi 



(c) , (c) / (z) . (z) 

Try m h uj h + m f Wf y mbw b + m f u; f 



(30) 



(31) 



according to Eqs. (fT8 |) - (f2T|) . 

Finally we obtain the quasi-one-dimensional effective 
Hamiltonian denoted by 



ff = ff b + ff f + ffbb + ffbf, 



(32) 



where these parts are obtained in Eqs. (|22j) . (|23j) . ([28 
and (j2"9"j) . 



B. Quantum- mechanical dynamics 



The renormalized coupling constants y b b and <?bf in 
Eqs. (|2"4"|) and (|2"5|) are determined as 



In quantum mechanics, dynamics is described as an 
initial value problem for the Schrddinger equation, 



ffbb 



47T 

1 



"Abb 



(26) 



i±\*(t))=H\9(t)) 



(33) 
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with the time-dependent state |\&(t)). Here the Hamil- 
tonian H is given in Eq. (1321) in the present paper. 

In order to obtain the state ^(t)) in Eq. (|33|). we need 
to give an initial state [^(to)) at the initial time to- In 
general, quantum-mechanical dynamics depends on the 
initial state. 

Although we can freely give the initial state in theory, 
it is seriously difficult to determine the initial state in 
actual experiments. That is because the state is in an 
extremely-huge configuration space in quantum many- 
body systems, and we can not obtain the whole infor- 
mation of the state experimentally. It is a fundamental 
problem in the many-body physics. 

In order to avoid this problem, we utilize the ground 
state |g) of the energy eigen-states, 



(H + H')\g)=E' g \g), 



for the initial state, i.e. 



|*(*o)> = |g>, 



(34) 



(35) 



and suddenly remove the additional potential part H' in 
Eq. ([34|) at t = to- In this scheme, the ground state 
can approximately be obtained in the experiments of the 
ultra-cold atomic gases owing to the cooling technique, 
and we can precisely control the additional potential by 
using the optical potential. 

In order to study the quantum-mechanical dynamics, 
we must observe, at least, a time-series of a physical value 
denoted by A(t) = (t)\A\i& (t)) for the corresponding 
physical value operator A. The physical value A(t) re- 
flects information of the time-dependent state fill ]. 

In the toroidal gases, two of the important physical 
values must be the toroidal densities n\,(0,t) and rif(0,t) 
defined as 

/>oo />oo 

n a (9,t) = dr c r c dr z <*(i)|^ p ) t (r)^ ) (r)|l'(i)) 

Jo J-oo 

(36) 

for a = b and f in the limited configuration space. These 
distributions directly reflect the diagonal elements of the 
one-body density matrices in the spatial representation 
and can experimentally be observed in the optical ab- 
sorption imaging method. 

The toroidal densities nb(9,t) and ri{(9,t) in Eq. (|36|) 
are reduced to 



Since the toroidal densities n,b(0, t) and ni(9, t) are pe- 
riodic real-valued functions for 9, they can be written as 
the Fourier series expansion, 

n a (0, t) = — - + E [a an {t) cos (nO) + b an (t) sin (n9)} 

Z7T ' — ' 

71=1 

(39) 

for a — b and f. Here the constant terms in Eq. ([59")) 
proportional to the particle numbers Nb and Nf of the 
bosons and fermions, respectively, because of the defini- 
tion of the particle numbers, 



2tt 

Y„ = / dO n a (6,t). 
'o 



(40) 



Here we assume that the particles mostly exist in the lim- 
ited configuration space, and dissipation to the residual 
space is neglected. 

In the present paper, we consider the even functions, 
b an {t) = 0, and focus on the lowest term, 

a al (t) = —d a (t), (41) 

where we introduce the deformation parameters 
1 [ 2 * 

d a (t) = Tr d9 n a (9,t) cos (42) 

M a JO 

for a — b and f. Here the parity of the toroidal densi- 
ties is determined by that of the initial state formed from 
the additional potentials in Eq. (|34|) because the original 
Hamiltonian H has the parity symmetry. In addition, 
the main excitation of the density fluctuation is also de- 
termined by the additional potentials because H also has 
the rotational symmetry. The lowest term must be the 
easiest mode to excite the toroidal gases. 

In order to excite such a deformation, we prepare the 
additional potential part H 1 in Eq. (|34|) as 

H'=J2 I dr ^P\r)V a cos 9^ a p \r) (43) 

a=b,f 

with the small potential strengths Vb and V{. It is reduced 
to 



H' = V h db + Vfdf 
with the deformation operators 



(44) 



»b(M)« E E 



k)$ 



K— — QO K, ' — — OO 



2tt 



and 

n f (M)« E E 



v v e i(K'-K)6 



K — — OQ K, — — OO 



2- 



-<*(*)|6+6^|*(t)) (37) 



(*(t)|4c«/|*(t)) (38) 



according to Eqs. ([TI ]) -(|T5 ]) . ((T5|) . and (HU). 



and 



1 



dh =2 ^ [ b l b ^+K-ih 

K — — OC 



di = ^ E ( c t c «-i 



«.t 



K — — OO 



(45) 



(46) 



according to Eqs. fTT ]) -([I3 ]) . ([15]). and HHJ). In fact, the 
additional potentials in Eq. (|4"51) deform the initial state 
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with the even parity, i.e., d a (to) 7^ and b an (t) = in 

E q . dnnD- 

Since 

d a (t) = (^(t)\d a \^(t)) (47) 

for a — b and f, the energy eigen- value problem for the 
initial state, [^(io)) = |g), in Eq. can also be in- 
terpreted as an energy minimum problem for H, not for 
(H + H'), under constrained conditions for db(to) and 
di(to) in the Lagrange multiplier method, where Vt> and 
Vf correspond to the Lagrange multipliers for tit, (to) an( i 
c£f(fo), respectively. 

By using the complete set of the energy eigen-states 
|n) for H, 

H\n)=E n \n) (48) 

and 

00 

EHH = 1 ' ( 49 ) 

n=0 

we obtain 

00 00 

=EE j {En ~ Em)t C«nm (50) 
n— m— 

with the coefficients 

Co. nm = (g\n)(n\d a \m)(m\g) (51) 

according to Eqs. (j3"3")l and (1551) . Thus we can calculate 
the time-development of d a (t) from the information of 
the energy eigen- values E n and eigen-states |n). Note 
that the deformation parameters d a (t) in Eq. (|50[) are 
real numbers, d*(i) = d a (t), because of the physical value 
condition, = d a and C* m „ = C anm . 

In the present paper, we predict the time-series of 
dh{t n ) and df(t n ) from the initial time to (= 0) to a final 
time t max as the regular interval data, i.e., 

t„ = nA (52) 

for the sampling number n = 0, 1, . . . , (N— 1) and interval 
time A = t max /(N— 1), and calculate the discrete Fourier 
transform, 

A*(o; m ) = A J2 e-^ l -d a {t n ) (53) 

n=0 

for a — b and f, w m = 2nm/(NA), and m = 
(-JV/2),(-JV/2+l),...,(JV/2). 
Note that the Fourier transform 

/oo 
dt e~ lult d a (t) (54) 
-OO 

is reduced to 

00 00 

D »H = E I! 2lC «»™*( £ » - £ ™ - W ) (55) 
n—0 m— 



according to Eq. (l50l) . In the contiguous limit, i.e., 
i max — > 00 and A — > 0, the discrete Fourier transform in 
Eq. (|53")l gets close to the Fourier transform in Eq. (|55l) . 

In general, how to visualize a large number of S func- 
tions, e.g., those in Eq. (I5"5j) . is a technical problem in 
quantum many-body physics. In fact, the selection of 
the visualization methods obviously affects detailed pro- 
files of the graphics of the calculational results. 

In the present paper, we select the method of the dis- 
crete Fourier transform in Eq. (|53[) because it can repro- 
duce that in actual experiments, where the observation 
data can be obtained as the time-series of db(t n ) and 
df(t n ). 

C. Calculational schemes 

Here we explain the calculational schemes used in the 
present paper. First we explain the numerical exact di- 
agonalization method as the main scheme in the present 
paper. Second we also explain the HF approximation 
used for comparison. 

1. Exact diagonalization method 

In the exact diagonalization method, we numerically 
solve the energy eigen- value problems in Eqs. (|34p and 
(l48l) by diagonalizing the Hamiltonian matrix in a lim- 
ited space in the matrix representation, where the diag- 
onalization is exactly computed except for numeric er- 
rors. The quantum-mechanical dynamics can be calcu- 
lated from the eigen- values and eigen-states according to 
Eqs. ([501) and (f5T ]) . 

As an advantage of this method, we can obtain ex- 
act solutions except for numeric errors and limitation of 
the matrix space. On the other hand, as a disadvantage 
of this method, we can not demonstrate the calculation 
for large-number and strongly-interacting systems, where 
the needed matrix space becomes extremely huge and 
overflows from our computer memories. 

In order to demonstrate the calculation, we decide the 
matrix space as follows. 

First we make the base vectors by using the Fock space 
for the boson and fermion operators b K and c K , where 
the particle numbers TVb and Nf are fixed. These base 
vectors work well except for the highly-deformed and/or 
strongly-interacting gases as a result. 

Second we limit the range of n as |«| < Kb max for the 
bosons and |k| < Kf max for the fermions. This limit can 
especially affect the interaction energies for the contact- 
type interactions in Eqs. (|28|) and (|29|) ; However, these 
effects on the dynamics are visually small except for the 
case of the strongly-interacting gases. 

Third we also limit the matrix space by eliminating 
the base vectors with the higher kinetic excitation ener- 
gies than a cut-off energy e c . Here the kinetic excitation 
energy of one of the base vectors means difference of the 
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expectation values of (£T b + iff) in Eqs. (|22|) and (f23|) 
between the base vector and ground state of the non- 
interacting gases. This limit is valid when the cut-off 
energy is much greater than the typical scales of the ex- 
citation and correlation energies, where the correlation 
energy is defined as difference between the interaction 
energy and mean-field energy. 



2. The Hartree-Fock approximation 

The HF mean-field approximation must be one of the 
most popular schemes to treat the ultra-cold atomic 
gases. Here, as a naming convention, the "Hartree-Fock" 
approximation is used for fermions, and the "Hartree" 
approximation is used for bosons. In the present paper, 
we use the "Hartree-Fock" approximation for the boson- 
fcrmion mixtures as the unified name. 

In the HF approximation, the configuration space of 
the state l'F(i)) is limited to that of the HF state denoted 
by 



|*hf(*)) 



iVi. 



&t(t) N t 

H4W|0) (56) 



with the vacuum |0) and redefined boson and fermion 
operators 



are neglected in this scheme. The validity of this scheme 
is discussed in the following sections. 

The single-particle wave-functions f K (t) and u nK (t) in 
Eqs. ([571 and (|58p are determined by the variational min- 
imization of the action 

Shf = J alt (# HF (t)|iA - ff|* HF (i)) (62) 

for the dynamics in Eq. (|33[) and of the energy 

E' UF = (VMt )\H + JT'|*HF(to)) (63) 

for the ground and initial state in Eq. Q34p with the nor- 
malization constraint in Eqs. (j59")l and (|60l) . 

After all, for the dynamics in Eq. (|33|) . we obtain the 
coupled time-dependent HF equations, 



j OO 

*—(*) = E 



dt 



K — — CXJ 



5bb N h — 1 (b)(HF) 
- - - -Pk'-k (*) 



9bf Jf)(HF),A 



and 

. dUr, 



dt 



'-it) 



E 



2m h R 2 2ttR N h 



JJW (b)(HF) m 



2m fj R 2 2?ri? 



(64) 

^tik' (^) 

(65) 



with the momentum distributions 



(57) 



(b) (HF) 



(t) = JV b 53 <^,(t)^_ K (i) 



(66) 



K — — OO 



and 



(58) 



K— — OO 



where the single-particle wave-functions ip K (t) and u nK (t) 
obey 



53 ^(tK(t) = i, 



U*n K (t)Un' K (t) = 5n 



and 



53 KMunK'it) = s K , K >. 



(59) 



(60) 



(61) 



n=l 



In the HF approximation, information of the many- 
body state is abbreviated into the single-particle wave- 
functions <p K (t) and u nK (t) in Eqs. (|57|) and (|55|) related 
to the one-body density matrices. As a result, the many- 
body correlations beyond the one-body density matrices 



and 

Nf oo 

P m( HF )(^53 53 < K ,(iK(K<- K) w. (67) 

n—1 k' — — oo 

Note that, if we replace the coefficient (A b — 1)/-/V b in 
Eq. with unity by neglecting , Eq. corre- 
sponds to the Gross-Pitaevskii equation. 

For the ground and initial state in Eq. (|34[) at t = to, we 
obtain the coupled HF equations for (p K (to) and u nK (to), 



E 



K?6 K K ' 



2m b i? 2 

fV — LXJ 

ffbb A b - 1 (b)(HF) 



Vb 



2tt.R JV b 



(*o) 



ffbf (f){HF) f , ^ 



y K (68) 



and 



E 



K — — OO 



2m fj R 2 



+ V £ 



5bf Jb)(HF) r . , 
+ 2^R P *'-* (to) 



(69) 



with the single-particle energies e^ b ^ and . In the HF 
ground state, according to Eq. (|56"|) . all of the bosons 
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occupy the lowest single-particle state, and the fermions 
exhibit the Fermi degeneracy. 

In the small amplitude excitations from the HF ground 
state, the HF approximation ag ree with the random 
phase approximation (RPA) [29l [37j described in ap- 
pendix [A] In the present framework, the amplitudes of 
the excitations are not limited in the small amplitude 
excitations, and we thus use the HF approximation. 

In addition, instability of the HF ground state can be 
verified by RPA because the instability can be seen in the 
small amplitude excitations as explained in appendix [A] 
The instability condition for the HF ground state is dis- 
cussed in the next section. 



III. STATIC PROPERTIES 

We here roughly explain static properties of the 
toroidal gases in the ground state as an important back- 
ground knowledge to study the dynamics. In Sec. IHI Al 
we explain effects of the quantum statistics and inter- 
particle interactions on the toroidal gases. In Sec. IIII Bl 
we also explain the static properties of the non-deformed 
toroidal gases without the deformation potentials and 
consider the instability condition of the HF ground state. 



A. Effects of quantum statistics and interactions 

In the non-interacting gases, all of the bosons occupy 
the lowest single-particle state, and the fermions exhibit 
the Fermi degeneracy, in common with the HF ground 
state. In fact, the ground state of the non-interacting 
gases agrees with the HF ground state. 

If we add the same deformation potential for the 
bosons and fermions to the non-interacting gases as de- 
scribed in Eq. (1341) , the bosons exhibit larger deformation 
than that of the fermions. That is because the fermions 
need more energy to deform the gas owing to the Fermi 
pressure. It is a typical difference between the bosons 
and fermions. 

Similarly the fermions are more insensitive to the inter- 
particle interactions than the bosons owing to the Pauli 
exclusion principle in two mechanisms as follows. First 
the fermion-fermion interaction is neglected in the model 
in Eq. (|3"2)) because of the Pauli blocking effect. Second 
the fermions has more kinetic energy than the bosons 
because of the Fermi degeneracy, and the motion of the 
fermions relatively reduces effects of the interactions. As 
a result, in a rough picture, the interactions chiefly affect 
the bosons, and the fermions reluctantly obey the bosons 
owing to the boson-fermion interaction. 

Although the three-dimensional gases can exhibit the 
collapse phenomena owing to the inter-atomic attrac- 
tions, the quasi-one-dimensional gases can not collapse in 
theory. That is because the quantum-mechanical kinetic 
energy (proportional to L~ 2 with the typical length size L 



of the gases for the small L) overcomes the attraction en- 
ergy (proportional to in the quasi-one-dimensional 
gases. Of course, it is a theoretical artifact derived 
from the limited configuration space for the quasi-one- 
dimensional model, and the gases can collapse in the 
full configuration space because of the three-dimensional 
attraction energy (proportional to L~ 3 ). The theoret- 
ical stability of the quasi-one-dimensional gases, how- 
ever, indicates realization of long-lived meta-stable states 
with the attractions in actual experiments. In fact, such 
states like the bright soliton are realized in the quasi-one- 
dimensional bosons with the attractions in the ultra-cold 
atomic gases. 

In the present paper, we thus theoretically consider not 
only the inter-atomic repulsions but also the attractions. 



B. The non-deformed systems 

In the non-deformed toroidal gases without the defor- 
mation potentials, the ground state must have the rota- 
tional symmetry in the finite quantum systems, where 
spontaneous symmetry breaking is forbidden. In fact, 
the ground state calculated in the exact diagonalization 
method has the rotational symmetry. 

Since the HF ground state is merely an approximation 
for the exact ground state, it may be unstable for the 
density fluctuation, and the instability of the HF ground 
state is usually related to symmetry breaking in general. 
Here it should be noted that the symmetry breaking for 
the mean-field theory is strictly different from the spon- 
taneous symmetry breaking in quantum field theory for 
the infinite matter systems. 

The instability condition of the non-deformed HF 
ground state is obtained from the lowest excitation en- 
ergy of the small amplitude excitations from the HF 
ground state, where the time-dependent HF approxima- 
tion is reduced to RPA as described in appendix [A] In 
fact, the excitation energy can go into an imaginary num- 
ber, where the corresponding fluctuation grows in time 
and indicates the instability. 

The instability appears only in the non-deformed sys- 
tems and disappears in the deformed systems with the 
deformation potentials. That is because the instability 
is linked to the symmetry breaking and the quasi-one- 
dimensional gases can not collapse in theory. 

The instability condition for the non-deformed systems 
can also be obtained from behaviors of the deformation 
parameters in the HF ground state in the deformed sys- 
tems with the small additional deformation potentials. 
Under the instability condition, the deformation parame- 
ters suddenly grow by addition of the infinitesimal poten- 
tials, and the gases exhibit the soliton-like profiles [34| . 
where the gases are spatially localized in order to balance 
the interaction and kinetic energies. 

In Fig.[TJ we show the instability condition as a typical 
demonstration, where the setting parameters and scales 
are same as those in Sec. II VI except for Vb = Vt = 0. 
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FIG. 1: The instability condition of the non-deformed HF 
ground state. The setting parameters and scales are same as 
those in Sec. IIVI except for Vb = Vi = 0. 



As shown in Fig. [TJ the boson-boson attraction can 
make the instability and induce the bright soliton for 
the bosons as a main mechanism of the instability. In 
particular, when g^t = 0, the instability occurs at 

gbb -' 2 TOb ^ b ~l) " g ^- (?0) 

This condition agrees with that in the time-dependent 
Gross-Pitaevskii equation, i.e., the Bogoliubov analysis, 
except for the difference term N^ 1 from the HF approx- 
imation. 

In Fig.[TJ the boson- fermion interaction also affects the 
instability and can induce the phase-separation-like pro- 
files for the repulsion, <?bf > 0, and the mixed soliton-likc 
profiles for the attraction, <?bf < 0, as the other mech- 
anism of the instability. As an aspect of this mecha- 
nism, the instability condition appears symmetrically for 
the boson-fermion repulsion and attraction as shown in 
Fig.ffl 

Although the instability of the HF ground state is a 
theoretical artifact in the HF approximation, it is deeply 
linked to the dynamics beyond the HF approximation as 
demonstrated and discussed in the following sections. 



IV. DYNAMICAL PROPERTIES 

In this section, we demonstrate the dynamics in some 
typical situations. In Sec. IIV A[ we first explain the set- 
ting parameters and scales needed to demonstrate the 
dynamics. In Sec. IIVB1 we show the calculational re- 
sults. 



A. Setting parameters and scales 

Since one of the purpose of the present work is to 
study effects of the quantum statistics on the dynam- 
ics, we select the symmetric parameters for the bosons 
and fermions, i.e., nib = rrif, Vb = V{, and N\, — Nf. In 
this case, both of the particles exhibit the same behaviors 
except for the effects of the quantum statistics. 

We also determine the scales as follows without loss 
of generality. In addition to ft = 1 set in the above sec- 
tions, we also select = mf = 1 for the mass scale and 
R/ (iVb + iVf) = I for the length scale. Note that we here 
select the inverse of the total density for the length scale, 
instead of R, because the total density is kept in the ther- 
modynamical limit, differently from R. In these scales, 
the typical absolute values of the coupling constants, |<? b b| 
and |<7bf|, in actual experiments become lO - ^ 3 ) depend- 
ing on detailed setup of experimental situations. 

We concretely set iV b = Nf = 5 and V h = V f = -0.001 
in the following demonstration. Here the precise value 
of Vb and Vf hardly affect the calculational results. The 
numbers TVb and Nf are chiefly determined by the limi- 
tation of our computer memories. In addition, it is theo- 
retically favorable for us to select an odd number for Nf. 
That is because the Fermi levels, which are degenerate 
double except for the ground state owing to the parity 
symmetry, are fully occupied for the odd number; Oth- 
erwise the single fermion on the Fermi levels exhibits a 
precise beat structure in the density distribution, and, 
after all, the calculational results for the odd number are 
clearer than those for the even number. 

In practice, the other parameters needed to demon- 
strate the dynamics are determined as follows. The cut- 
off momenta are determined as K b max = 4 for the bosons 
and Kf m ax = 5 for the fermions, where the imbalance 
for the bosons and fermions is due to the Fermi degen- 
eracy only for the fermions. The cut-off energy is de- 
termined as e c = 0.14. Here the cut-off energy is much 
greater than the typical scales of the effective excitation 
energies, which are shown in the results in the following 
figures, and absolute values of the correlation energies, 
which are comparable with or smaller than that of the 

mean-field energy, |e m . f .| - \g$\Ng/(4<irR) « 0.0078, in 
the weakly-interacting gases studied below. 

In addition, we also determine the sampling number n 
and maximum time i ma x in order to obtain the discrete 
Fourier transforms £>b and Df in Eq. as N = 2 11 
and i max = 20 x N. Here it should be noted that de- 
tailed profiles of the discrete Fourier transforms depend 
on the sampling parameters, N and i max , and these pa- 
rameters may have some practical limitations in experi- 
ments, e.g., the life-time of the gases, thermal damping, 
and so on. Thus the detailed profiles may not be serious 
in the present paper, and we here focus on the qualita- 
tive profiles robust over a large region of the sampling 
parameters. 
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B. Calculational results 

Here we show the calculational results when (fof < 
in Sec. IIVB1I and when g hi > in Sec. IIVB21 Some 
theoretical interpretations of the results are discussed in 
the next section. 

In the present paper, we especially focus on the boson- 
boson attraction case, (fob < 0. That is because the 
boson-boson attraction induces a qualitative change of 

(c) 

the dynamics around (fob = ffbb related to the instability 
of the HF ground state in Fig. Q] as demonstrated below. 



1. The boson-fermion attraction case 

First we show the results for the (fob = case in Fig. [2] 
in order to demonstrate effects of the boson-fermion at- 
traction on the dynamics. Second we show the results 
for the (fof = case in Fig. [3] in order to demonstrate ef- 
fects of the boson-boson interaction. Third we show the 
results for the (fob < and <fof < case in Fig. H) 
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FIG. 2: (Color online) The absolute values of the discrete 
Fourier transforms, \D^\ [(a) (c)] and \D{\ [(b) (d)], when 
gbb = and flttf = [(a) (b)] and -\g$\ [(c) (d)]. The 
top and bottom columns represent the results in the exact 
diagonalization method and HF approximation, respectively. 



In Fig. [2 we show the absolute values of the discrete 
Fourier transforms, |-D b | and \Df \ in Eq. (1531) . when (fob = 
0. Note that the Fourier transforms are even functions, 
Db(-u> m ) = D h (uj m ) and D { (~uj m ) = D f (w m ), because 
of the real- valued d b (t) and df(t), and visually zero in the 
upper outside of the plotted range in Fig. [5] and also in 
Figs. EE 

In Figs.^a) and(2fb), the results in (fob = (fof = ex- 
hibit single sharp peaks at co 0.0050 for the bosons 
and co ~ 0.0250 for the fermions, respectively. Here 



the results in the exact diagonalization method and HF 
approximation plotted in the top and bottom columns, 
respectively, are visually same except for the negligi- 
ble numeric errors. The positions of the peaks agree 
with the first excitation energies linked to the quantum 
statistics, i.e., l/(2mb-R 2 ) = 0.005 for the bosons and 
(3 2 - 2 2 )/(2™ f i? 2 ) = 0.025 for the fermions. The differ- 
ence of the vertical scales for the bosons and fermions 
indicates that of the oscillation amplitudes determined 
by the deformation parameters in the initial state, i.e., 
the effect of the Fermi pressure explained in Sec. IIII Al 

In Figs.^c) and[^d), the results in <fof = — | <?b»t? I ex- 
hibit almost the same structure in Figs. [2ja) and EJb) 
except for a new peak only for the fermions. The new 
peak is exactly analogous to the peak for the bosons in 
the shape and position. According to the HF approxima- 
tion in Eqs. and (f6"5"j) . the new peak for the fermions 
can be induced by the bosons as a forced oscillation ow- 
ing to the boson-fermion attraction. On the other hand, 
the effect of the fermions on the bosons is visually small. 

(b)g„= -0.5|g tl "| 
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FIG. 3: (Color online) The absolute values of the discrete 
Fourier transforms, \Db\, when pbf = and <?bb = (l/2)]</^j| 
(a), (-1/2)1^1 (b), -\g<£\ (c), and (-3/2)|«$| (d). The 
top and bottom columns represent the results in the exact 
diagonalization method and HF approximation, respectively. 



when (fof = 0, where \D{ 



In Fig. |31 we show |D b | 
same as that in Fig. 

In Figs. OUa) and[3]Jb), the results in gbb = (1/2) j^^b I 

and (— l/2)|(7kb | exhibit single peaks at to « 0.0061 and 
0.0036, respectively. Here the results in the HF approx- 
imation agrees with those in the exact diagonalization 
method because of the weak interactions. The heights 
of the peaks in Figs. Eta), [2fa), and E{b) reflect those 
of the oscillation amplitudes according to the Parseval's 
theorem. 

In Fig. EIc), the results at the instability boundary of 

(c) 

the HF ground state, gbb = <fob ' > exhibit a major peak at 
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u « 0.0021 and minor peaks near u « 0.0037. The HF 
approximation does not reproduce the minor peaks. 

In Fig. EJd), the results in gbb = (— 3/2)|g^| ex- 
hibit a major peak at w « 0.0012 and a minor peak at 
u) R3 0.0032. Here the results in the HF approximation 
exhibit quite different behaviors from the exact results 
and indicate the dynamics of the soliton. In fact, the re- 
sults in the HF approximation exhibit a dominant peak 
at u = 0, corresponding to a stationary soliton, and a 
minor peak at u) ~ 0.0066, corresponding to surface os- 
cillations of the soliton, where the minor peak is also 
slightly seen in Fig. [3fc) at ui 0.0063. 
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FIG. 4: (Color online) The absolute values of the discrete 
Fourier transforms, |D b | [(a) (c)] and |Df| [(b) (d)], when 

9bl = -|ffbb| and 9bb = -\g b b\ K a ) ( b )l and (— 3 / 2 ) Iffbb 1 1 
[(c) (d)]. The top and bottom columns represent the results 
in the exact diagonalization method and HF approximation, 
respectively. 

In Fig. 01 we show |Z>b| and \D{\ when = — |<? b b| 
and gbb < 0. These situations correspond to those in 
Figs. H^b), G3c), and [2Id) except for the boson- fermion 
attraction. 

I (c) I 

Here both of the results in <?bb = — fl^b anc ^ 
(—3/2) | I exhibit almost the same structure in 
Figs. HJb), [3jc), and[3][d) except for new peaks only for 
the fermions in analogy with the case in Fig. [3J The 
mechanism of the forced oscillation must not be affected 
by the boson-boson interaction as demonstrated here. 



2. The boson-fermion repulsion case 

Here we show the results when g^f > and comment 
on symmetric properties for the boson-fermion attraction 
and repulsion. 
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FIG. 5: (Color online) The absolute values of the discrete 
Fourier transforms, Dbj (a) and \Dt\ (b), when t/bf = |ffbt>| 

I (c) I 

and #bb = — <?j,b i " ^ and bottom columns represent 

the results in the exact diagonalization method and HF ap- 
proximation, respectively. 



In Fig. El we show \D^\ and |£>f| when g^i — 



I ( c )| 

Kb I 

I (c) I 

and gbb = — \ .<?bb I • This situation agrees with that in 
Figs. IDJa) and SJb) except for the sign of <?bf • 

As shown in Figs. |4|a), and|4]T>), |A>| an d |A| are 
mostly independent of the sign of gb{ in a wide region 
of the parameters. The symmetry for the boson-fermion 
attraction and repulsion is also shown in the instability 
condition of the HF ground state in Fig. [T] 

In the interpretation of the forced oscillation explained 
above, the sign of the external force merely affects the 
phase of the oscillation. In addition, the effect of the 
sign of gbf on the initial state merely affects the phase 
of the initial position of the oscillator. As a result, the 
symmetry can appear in the results. This interpretation 
can also be confirmed by viewing of the time series of 
df(t) before the Fourier transform. 



V. DISCUSSION 

In this section, we discuss the results in Sec. IIVI in 
some theoretical points of view. In Sec. IV A| we discuss 
the effects of the quantum statistics and interactions on 
the dynamics by viewing of the excitation energy spectra. 
In Sec. lVBl we also discuss the correlation effects beyond 
the HF approximation. 



A. Excitation energy spectra 

According to Eqs. ([50]) and (|55|) . the excitation energy 
spectra, E n — E = fl n obtained in Eq. f]48|) , are directly 
linked to the dynamics. 

In particular, the small amplitude excitations from the 
ground state are mostly construct from the n = or 
m = parts in Eqs. ([50]) and (|55|) . where these parts have 
the excitation energies 17 m or f£„, respectively. That is 
because the absolute values of the coefficients C a o m and 
Cano m Eq. (|5ip are much greater than those of C anm 
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for n ^ and to 7^ in the small amplitude excitations, 
i.e., |g) ~ |0). 

In the present scheme, the amplitudes of the oscilla- 
tions are determined by the interactions, gbb and <?bf, and 
deformation potentials, Vb and Vf, in the initial state. 

In particular, for the small deformation potentials 
demonstrated in the previous section, the amplitude for 
the fermions is small, and that for bosons is chiefly deter- 
mined by the boson-boson interaction, gbb, because the 
boson-fermion interaction has small contributions to the 
deformation owing to the Fermi pressure as explained 
in Sec. IIIII In a rough estimation, the oscillations for 
bosons can be treated as the small amplitude excitations 

(c) 

when gbb > g^b corresponding to the stability of the HF 
ground state in Eq. (|70]) . 

Here we first discuss the gbf = case, and then com- 
ment on the gbi ^ case. 



0.03 




0.005 



the bosonic spectra. Note that each of the spectra is de- 
generate, at least, double owing to the parity symmetry. 

Here it should be noted that the quantum statistics 
directly affects the excitation energy spectra linked to 
the dynamics. The results in RPA (denoted by the solid 
line) are also plotted in Fig. [5J The left boundary of 
the spectra in RPA (denoted by the thin-solid separa- 
tor line) indicates the instability boundary of the HF 

(c) 

ground state, i.e., gbb = g^b ■ Here the results in RPA 
reproduce the exact ones for the lowest fermionic spec- 
tra exactly and for the lowest bosonic spectra except for 
the neighborhood of the instability boundary. Note that, 
in the particle-hole RPA, we can not obtain the multi- 
particle-multi-hole spectra, as shown in the bosonic spec- 
tra in Fig. [5J because of the ansatz for the Q operators 
in Eq. (|A%|) . 

In the small amplitude excitations, the dominant peaks 
appear at the lowest spectra as explained above. It is 
confirmed by the results in Figs. HJa), Hfb), [21a), and 
Mh). 

Moreover we can find out that, in the non-small ampli- 
tude excitations in Figs.[3][c) and[3]Jd), the major peaks 
also appear at the lowest spectra, and the minor peaks 
correspond to the mixing of the lowest and next-lowest 
spectra. In addition, we also see that the invalidity of the 
HF approximation for the dynamics is directly linked to 
the instability of the HF ground state as shown in the 
lowest spectra in Fig. [5J 

Lastly we comment on the <?bf 7^ case. 

The boson-fermion interaction slightly mixes the exci- 
tation energy spectra in Fig. [5J In particular, the mixing 
of the lowest bosonic and fermionic spectra is perturba- 
tive because of the energy difference of the lowest bosonic 
and fermionic spectra derived from the quantum statis- 
tics. As a result, this perturbation produces the forced 
oscillation as shown in Figs. HJc), HJd), [H and [5] 
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9bb 



FIG. 6: (Color online) The excitation energy spectra, Q n , 
when gbi = 0. The dashed and solid lines indicate those in 
the exact diagonalization method and RPA, respectively. The 
left boundary of the spectra in RPA denoted by the thin-solid 
separator line indicates the instability boundary of the HF 
ground state. 



B. Correlation effects 

Here we discuss the correlation effects on the dynamics. 
We first introduce indicator parameters for the many- 
body correlations and the HF condition used in many- 
body physics, and then discuss the correlation effects by 
viewing of the calculational results of the indicator pa- 
rameters. 



1. The Hartree-Fock condition 



In Fig. [6j we show the excitation energy spectra, fl n , 
in the gbt = case, where the bosons and fermions are in- 
dependently moving as shown in Figs. [31 [3Ja), and[2Jb), 
and the spectra (denoted by the dashed line) are com- 
posed by the bosonic and fermionic spectra independent 
of each other. In fact, we can find the lowest fermionic 
spectra independent of gbb with D, = 0.025 in Fig. [51 and 
the other spectra below the lowest fermionic spectra are 



The indicator parameters are defined as 



for bosons and 



Ob 



CTf 



fr[/>f] 
N { 



(71) 



(72) 
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for fermions with the one-body density matrices 

( Ph ) KK ,=(g\btb K ,\g) (73) 

and 

(«)™< = (g|4<v|g) (74) 

for the ground state |g). Here the indicator parameters 
are independent of selection of the base set for the density 
matrices because the traces in Eqs. (|7T|) and (|72[) are 
independent of it. In addition, < Ob < 1 and < erf < 
1 because of the definitions in Eqs. (|7T|) and (fT2]) and 
positive-definite density matrices. 

In the HF approximation, we can obtain cr^ = Of = 1, 
which is called the HF condition. That is because the 
eigen-values of the density matrices agree with the occu- 
pation numbers of the single particles in the HF approx- 
imation. 

In general, the values of at, and at are less than unity 
owing to the many-body correlations beyond the HF ap- 
proximation. Thus they are usually used for the indica- 
tors of the correlations in many-body physics. 

2. Results 



In Fig. we show ah and <7f . 




9bb 9bb 

FIG. 7: (Color online) The indicator parameters Ob (a) and <7f 
(b). The solid and long-dashed lines indicate the results when 
Vh = Vf = and (?bf = and — \g hh | , respectively. The dashed 
and dotted lines indicate the results when Vh = Vf — —0.001 
and gbf = and — <? bb , respectively. The thin-solid separator 

lines indicate <7bb = S b b • 

First we focus on the results when Vb = Vf = and 
<7bf = (denoted by the solid lines). In this case, e>f = 1 
owing to the non-interacting fermions. On the other 
hand, <7b decreases from unity as |<7bt>| increases. In par- 
ticular, (7b noticeably decreases around <?bb = fl^b (^ e ~ 
noted by the thin-solid separator lines). 

Second we focus on the results when ghf = — | (de- 
noted by the long-dashed and dotted lines). The boson- 
fcrmion interaction slightly decrease the values of at, and 
o-f. As shown in Fig. the correlation effects of the 
boson-fermion interaction are much smaller than those 
of the boson-boson interaction. 

Third we finally focus on the results when Vh — Vf = 
—0.001 (denoted by the dashed and dotted lines). The 



small deformation potentials greatly increase the values 
of ah and at around the instability boundary. 

In general, the many-body correlations increase as 
|<7bb| an d \9bt\ increase as shown in Fig. [TJ However, their 
effects on the dynamics are visually small, i.e., the HF 
approximation is valid, except for the instability region 
as demonstrated in Sec. HVBl 

In the instability region, the many-body correlations 
must be related to the rotational symmetry. That is be- 
cause, as shown in Fig. [71 the small deformation poten- 
tials greatly reduce the correlations by breaking the ro- 
tational symmetry in the instability region; In the other 
region, the deformation potentials hardly affect the val- 
ues of (7b and <7f. 



3. Discussion 

Here we discuss three points about the effects of the 
many-body correlations on the dynamics. 

First we discuss a relationship between the many-body 
correlations and rotational symmetry in the instability 
region in terms of the HF theory. The relationship is 
demonstrated in Fig. [7] as explained above. 

In the HF theory, the instability induces the deformed 
HF ground state with an infinitesimal deformation po- 
tential called a seed of the deformation. It is consistent 
with the great reduction of the many-body correlations 
with the small deformation potentials in the instability 
region in Fig. [7] 

The deformed HF ground state can not directly be used 
for an approximation of the exact ground state. That is 
because the deformation direction of the deformed state 
is determined by the artificial seed, and the different 
seed produces the different state. Furthermore all of the 
produced states are degenerate because of the rotational 
symmetry. In principle, the ground state must have the 
rotational symmetry in the finite quantum systems. 

As an extended HF theory, we can consider a superpo- 
sition state of all of the deformed HF ground states and 
determine the coefficients by estimation of the full energy 
beyond the mean-field energy. This scheme is called the 
generator coordinate method or projection method and 
studied well in nuclear physics for deformed nuclei [29 ], 

In the projection method, we can obtain the approx- 
imate states for the ground state and some low-energy 
excited states, and the rotational symmetry of the states 
is restored. 

In this sense, the many-body correlations in the insta- 
bility region can be interpreted as the superposition of 
the deformed HF ground state for the restoration of the 
broken symmetry. In addition, the bosonic low-frequency 
modes in the dynamics in Sec. IIVBI can also be inter- 
preted as the restoration modes of the broken symmetry 
in terms of the projection method. 

Second we discuss another relationship between the 
many-body correlations and macroscopic behaviors of the 
toroidal gases in the thermodynamical limit, where the 
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macroscopic behaviors can be described in quantum field 
theory except for the finite-size effects. 

According to quantum field theory, the ground state 
can exhibit the spontaneous symmetry breaking with the 
inequivalent vacua and Nambu-Goldstone mode in the 
thermodynamical limit [38j . where quantum field theory 
must be exact, unlike the HF theory. 

The symmetry broken states exhibit macroscopic ob- 
jects and can not transfer to each other by any quantum- 
mechanical operators. In fact, the dynamics in Sec. II V 51 
must get close to the macroscopic behaviors, where the 
dominant oscillation mode is almost at u) — owing to 
the correlation effects and indicates the immobile object 
in a finite time 42]. In other words, the dynamics must 
be dominated by the finite-size effect. 

In terms of quantum field theory, the many-body cor- 
relations related to the symmetry can be interpreted as a 
precursor of the spontaneous symmetry breaking, or the 
phase transition, in the finite quantum systems. 

Third we discuss a mechanism of the many-body cor- 
relations in quantum mechanics for the finite systems. 

In quantum mechanics, the many-body correlations 
related to the symmetry must be due to existence of 
quasi-degenerate states. That is because the perturba- 
tive deformation potentials greatly affect the correlations 
as shown in Fig. [Jj and the perturbation can not induce 
such a great effect on the states except for the quasi- 
degenerate states according to the perturbation theory. 

Since a degeneracy in quantum mechanics must be 
based on a symmetry and its irreducible orthogonal 
states, the quasi-degeneracy must also be based on the ro- 
tational symmetry and approximately-orthogonal states, 
(g(#d) g(0d')) ~ Sg d _g d ', where we introduce the ground 
state, g(#d)) for H + H'(9d) with the deformation direc- 
tion #d (= for the initial state in the dynamics). The 
approximate orthogonality can clearly appear near the 
the thermodynamical limit corresponding to the phase 
transition in quantum field theory and also gradually ap- 
pear in the small-number systems. 

In the finite quantum many-body physics, the quasi- 
degeneracy must be important for the correlation effects. 

Although the above three points are based on the dif- 
ferent theories in mathematics, the findings must be con- 
sistent in terms of physics. Thus the studies on the dy- 
namics in the cold atoms must contribute to those on the 
theories. 



In Sec. IIII1 we explain the static properties of the 
toroidal gases as a background knowledge. 



In Sec. IIV1 we demonstrate the dynamics in the sym- 
metric situations for the bosons and fermions. In partic- 
ular, we show the quite different behaviors of the bosons 
and fermions owing to the quantum statistics, and also 
demonstrate the weak correlation effects in the stabil- 
ity region and strong correlation effects with the low- 
frequency modes in the instability region. In addition, we 
reveal the strong influence of the bosons on the fermions 
as a forced oscillator. 

In Sec. El we discuss the effects of the quantum statis- 
tics and many-body correlations on the dynamics by 
viewing of the excitation energy spectra and indicator 
parameters for the correlations. 

The present study must be linked to actual experi- 
mental situations; However, although the toroidal bosons 
have been realized and studied well, the experimental 
studies on the boson-fcrmion mixtures of the toroidal 
gases have not been demonstrated yet. We hope real- 
ization of the studies in future experiments. 

Our scheme in the present paper must be applicable 
to various studies on the dynamics and many-body cor- 
relations in the cold atoms. In particular, studies on 
the many-body correlations may be important in the 
recently-realized systems with the dipolar interaction (39j 
and spin-orbit interaction [|o| because these systems may 
have rich phase structures related to the correlations. 

Lastly we comment on the previous studies on the col- 
lective excitations of the boson- fermion mixtures [3, . 
These studies cover the three-dimensional gases and are 
based on the HF approximation. Although the present 
study covers the quasi-one-dimensional gases and is be- 
yond the HF approximation, the results in these studies 
show some analogous behaviors, e.g., the different be- 
haviors of the bosons and fermions, the influence of the 
bosons on the fermions as a forced oscillator, and so on. 
These behaviors may be independent of the dimension- 
alities and many-body correlations. 

In future works, we should study the asymmetric situ- 
ations for the bosons and fermions, which are not demon- 
strated in the present paper and must be treated in the 
same formulation. In addition, we should also study the 
toroidal gases with the very strong interactions, which 
are not treated in the present paper. 



VI. SUMMARY AND PERSPECTIVE 

In the present paper, we study the dynamics of 
the quasi-one-dimensional boson-fermion mixtures of the 
toroidal gases in the exact diagonalization method. Es- 
pecially we focus on the effects of the quantum statistics 
and many-body correlations. 

In Sec. mi we formulate the model and dynamics, and 
also explain the calculational schemes. 
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Appendix A: Random phase approximation 



and 



In RPA [29j], we approximately obtain a part of the 
excitation energies 



E n — E 



and Q operators 



Qn = |0)<n| 



(Al) 



(A2) 



for the quantum number n (> 1) in Eqs. (j48|) and (j49|) 
by using the variational method in a variational ansatz 
for the Q operators to minimize the excitation energies. 

Because of the orthonormal relations, the Q operators 
must satisfy 



(o\Q„Ql\o) = (o\[Q n ,Ql]_\o) = s n , 



(A3) 



where the commutation relation in the middle part is a 
convention. Note that the Q operators are not the boson 
operators owing to the definition in Eq. (|A2|) . 

The Q operators determine the excitation energies in 
Eq. (|aT1) as 



(0 


Qn 


[h,qi 


]_ 


0) 


(0|Q„qJ, 


0) 



(0\[Qn,Qt]_\0) 



where the right part is also a convention. 

In the particle-hole RPA, we assume the Q operators 

as 



Xfttbo - Y k b\b k 



E^ 

1 v 

E E (< h 4 



pC-h Vn C'h 



|p|<Kf \h\>K { 



W) (A5) 



with the Fermi momentum Kf and variational parameters 



k vk X p'\ and y% h . 



In addition, the ground state |0) in Eqs. (|A3|) and (IA4|) 
is approximated by the HF ground state |0hf), where 

&A;|0hf) = c p |0hf) = c 1|0hf) = and so on. Then the 
particle-hole RPA agrees with the time-dependent HF 
approximation in the small amplitude excitations from 
the HF ground state. That is because the HF state in 
Eq. (|56[) in the small amplitude excitations can be con- 
structed from the particle- hole fluctuations in Eq. (|A5I) . 

As a result, we obtain the variational equations, 5fl n — 
0, denoted by 



with 



Az n ^l n z Ti 



{X k ; k £ 0} 
{Y n k : fc^O} 

Xf] H < Kf, \h\ > Kf 

yP h ; \p\<Ki, \h\>Kt 



(A6) 



(A7) 



tkp'h' 
A U 

*kp' h' 

I I ) I — J I -' J I 24 

c ' ' ' '' c ' 1 ° ' ' ' ' ' ^ r Aphp'h' 

^34 

A php'h' 
^44 

' (A8) 

where the superscripts in Eq. (|A8|) vary in the same 
ranges in Eq. (1A7[) in the same order, and fc, p, and h 
indicate the column components; fc', p', and h! indicate 
the row components. The Hermitian conjugate becomes 




(A9) 



with 



D 



{S kM } {0} {0} {0} 

{0} {-6 k , k ,} {0} {0} 
{0} {0} {S P y6 hJl ,} {0} 
{0} {0} {0} {-5 p ,p,5 hM ,} . 

(A10) 

because A* = DAD, D 2 — 1, and fl* n = tt n . 

The matrix elements in Eq. (|A8I) are obtained as 



a kk 

/i n — 

\ kk' 

Akp'h _ 

Akp'h' _ 
^14 

tphk' _ 
-^31 — 

Aphk 

A 32 — 

Aphp'h' 
^33 



-A, 



h:h: 



+ 9bh{N h - 1) 



2m h R 2 
= gbb (Nb - 1) <5 ,(fc+fc')' 

-A kph = gbi\/Nb$p>,(h>+k), 

-A kph = gbf\/~Nb5h',(p'+k)-, 

-A P ^ k — gbf\fNb&p,(h+k'), 

-A^ = gbf\fNh&h,(v+k'), 
?-h 2 



APhp'h' 



2m f R 2 



&p,p'3h,h' ; 



j^php'h' _ j\P h P h ' — Q 



43 



(All) 



The orthonormal relations in Eq. (|A3[) become 

4 Dz nl = 5 n , n , (A12) 

and can be kept in Eq. (|A6|) except for degeneracy be- 
cause of the real- valued energy Q„ and Eq. (|A9|) . 

Eq. (|A6|) has pairs of the symmetric solutions for the 
positive and negative energies (denoted by 0„ = — fi_ n > 

0) owing to the time-reversal symmetry. Of course, the 
physical solutions must be the positive solutions (for n > 

1) because of the definition in Eq. (|A1[) . 

Although the physical lowest eigen- value fli in 
Eq. (IA4|) must be a real positive number for the exact 
ground state |0), it can go into a pure imaginary number 
for the HF ground state |0hf)- The imaginary solution 
indicates the instability of the HF ground state for the 
particle-hole fluctuations. 
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